double dist( double x1, double y1, double x2, double y2 )
{
    double s1,s2;
    s1=pow((x1-x2),2);
    s2=pow((y1-y2),2);
    return sqrt((s1+s2));
}